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Searching for the signature of the violation of chiral charge conservation in solids has inspired 
a growing passion on the magneto-transport in topological semimetals. One of the open questions 
is how the conductivity depends on magnetic fields in a semimetal phase when the Fermi energy 
crosses the Weyl nodes. Here, we study both the longitudinal and transverse magnetoconductivity of 
a topological Weyl semimetal near the Weyl nodes with the help of a two-node model that includes 
all the topological semimetal properties. In the semimetal phase, the Fermi energy crosses only 
the 0th Landau bands in magnetic fields. For a finite potential range of impurities, it is found that 
both the longitudinal and transverse magnetoconductivity are positive and linear at the Weyl nodes, 
leading to an anisotropic and negative magnetoresistivity. The longitudinal magnetoconductivity 
depends on the potential range of impurities. The longitudinal conductivity remains finite at zero 
field, even though the density of states vanishes at the Weyl nodes. This work establishes a relation 
between the linear magnetoconductivity and the intrinsic topological Weyl semimetal phase. 


I. INTRODUCTION 

Searching for the violation of chiral charge conservation in solids started with Nielsen and Ninomiya’s proposal 
in 1983 [1], in which the chiral charge is not conserved in a ID system of two bands with opposite chirality. To 
simulate the ID chiral bands, they proposed to use the lowest Landau bands of a 3D semimetal, and expected that 
the longitudinal magnetoconductance becomes extremely strong. Recently, thanks to the discovery of a number of 
realistic materials of topological semimetals [2-21], there is a growing passion on their electronic transport [22-30] 
and signatures of the chiral anomaly [31-35]. 

Earlier theories on the longitudinal magneto conductivity arrived at various results [36-43]. In the semiclassical 
limit, where the Landau levels are not well formed, a positive magnetoconductivity was predicted [36, 37], and is 
under intensive experimental investigation recently [44-51]. In the semiclassical approaches, the Fermi energy should 
overwhelm the relaxation rate, not exactly at the Weyl nodes. The magnetoconductivity is also obtained by 
modeling the disorder as long-range charged impurities in the quantum limit [41]. In the scenario similar to that 
proposed by Nielsen and Ninomiya, different results have been obtained so far, depending on models and treatments 
[37-41]. Literally, a semimetal must have a Fermi energy crossing the Weyl nodes. Nevertheless, little attention is 
paid to the magneto conduct ion in the exact semimetal phase. More importantly, the Weyl nodes always appear in 
pairs. The intrinsic connection of the Weyl nodes and the inter-node scattering are two factors to affect the transport 
properties because the chiral anomaly occurs between two Weyl nodes. 

In this work, we start with a two-node model to investigate both the longitudinal and transverse magnetocon¬ 
ductivity of a Weyl semimetal near the Weyl nodes. The model describes a pair of Weyl nodes, and is solvable in 
the presence of magnetic fields. The scattering potential of the impurities is modeled by using a random Gaussian 
potential, in which the range of potential may vary in realistic materials. As long as the potential range is finite, 
we show that the longitudinal magneto conductivity is positive and linear in magnetic field, giving rises to a negative 
magnetoresistance. As the held goes to zero, we have a hnite minimum conductivity, even though the density of 
states (DOS) vanishes at the Weyl nodes. In the transverse magnetoconductivity, we hnd a crossover from linear B 
dependence in the short-range potential limit to 1/B dependence in the long-range potential limit. 

The paper is organized as follows. We hrst introduce the two-node model and present the solution of the Fermi 
arcs in Sec. 11. We then discuss the magnetic held-induced DOS and chiral anomaly in Sec. III. In Secs. IV and V, the 
formulas and analysis of the longitudinal and transverse magnetoconductivities in the presence of the random Gaussian 
potential are given, respectively. We also discuss the paramagnetic Weyl semimetal in Sec. VI. The conclusions are 
given in Sec.VII. The detailed calculations are provided in Appndix A-G. 
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FIG. 1: The energy spectrum of the bulk states (a) and the surface states at the y=Q surface (b) of the topological Weyl 
semimetal, (c) Both the bulk and surface states, (d) The real-space schematic of the topological Weyl semimetal and its 
surface states. We asusme that the x and z directions are infinitely long. The lines with arrows in the i/=0 (red) plane indicate 
that the chiral surface states travel along only one direction. In contrast, a topological Dirac semimetal hosts helical surface 
states that travel along both the x and —x directions. Parameters: = 0.1/nm, M = 5 eV-nm^, and A = 1 eV-nm. 


II. TWO-NODE MODEL OF WEYL SEMIMETAL AND FERMI SURFACES 


We describe the topological semimetal with a two-node model [40], 

H{k) = A{kx(Tx + kyGy) + M{kl - kl - k^ - kl)a^, (1) 

where (Jx,y,z are the Pauli matrices, k = (kx, ky,kz) is the wave vector, and A, M and k^ are model parameters. The 
model is in the topological semimetal phase and describes a pair of 3D gapless Dirac cones, with two Weyl nodes 
located at k = (0, 0,±fcc) in momentum space [see Fig. 1(a)]. The topological properties of the two-node model can 
be examined by the Berry curvature, Chern number and Fermi arcs [40, 52]. A topological semimetal has the kz- 
dependent topologically protected surface states for a specific kz between the two Weyl nodes. This is demonstrated 
by a non-zero Chern number as a function of /cz, Nc{kz) = sgn(M) for \kz\ < kc, and 0 for \kz\ > kc- According 
to the bulk-edge correspondence, there exist surface (or edge) states around the surfaces parallel to the z direction. 
The solution of the surface states can be found from the two-node model explicitly by following the solution to the 
two-dimensional modified Dirac equation that describes the quantized anomalous/spin Hall effects [53]. 

Suppose we have a semi-infinite system in the half plane y ^ 0 with open boundary conditions and with translational 
symmetry along the x and z directions, as shown by Fig. 1(d). kx and kz are still good quantum numbers but ky is 
replaced ky = —idy in the Hamiltonian (1). We can assume a trial wavefunction for each set of kx, kz as 


V'A = e 


ikxX-\-ikzZ 



( 2 ) 


Substituting the trial wavefunction into the eigen equation H{kx,—idy,kz)ip\ = Eip\, we have the secular equation 
for the eigen energies 


det \H{kx, -iX, kz) - E\ = 0, 


(3) 
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which gives four solutions of X{E), denoted as /3 Aq with /3 = ±, a = 1, 2, and 

^a[^) — I 2M2 + 2M2 

where Ak = — k^. Each f3Xa{E) corresponds to a spinor state 

i’aP = 


M(Afc + Xl) + E 
A{kx + /3Aq) 


or 


A{kx — /3Xa) 
-M{Ak + Xl)+E 


ijjalS = 

The general wavefunction in the y direction can then be written as a superposition of the spinor states 


^k^,kAE,y)= ^ , 

a=l,2,/3=± 


( 4 ) 

( 5 ) 

( 6 ) 

( 7 ) 


where the eigen-energy E as well the coefficients Cap are to be find from the boundary conditions. 

We now apply the open boundary conditions: (E, —oo) = (7?; 0) = 0. The former condition 

'^k^,k^(E,—oo) = 0 requires that t]/ contains only the terms with positive /? and Re(AQ) > 0 (or negative /3 and 
Re(Aa) < 0) and thus Ci- = C 2 - = 0. According to Eq. (4), this condition can only be satished in two cases: (i) 
A 1.2 > 0, and (ii) Ai _2 = ib with a,b > 0 (Note that Ai = A 2 corresponds to a trivial case). The later condition 
'^k^,k^{E,0) — 0 then gives a secular equation 

det|V'i-H 'ip2+\=0 (8) 

to determine Ci+ and C 2 +- Substituting Eqs. (5) and (6) into Eq. (8), respectively, and considering Ai X 2 , we 
arrive at 


E— —M{Ak — X 1 X 2 ) + Mkx{Xi + X 2 ), (9) 

E = Af(Afe — A 1 A 2 ) -l- Mkx{Xi -\- A 2 ), (10) 

which lead to 

A 1 A 2 = Afc. (11) 

On the other hand, according to Eq. (4), we have 

A?A2 = {A^kl- E^)IM^ + Al 
Xl + X^ = A^/M^-2Ak. 

Substituting Eq. (11) into Eq. (12), we have immediately 

E^ = A^kl- 

Using Eqs. (11) and (13) and keeping in mind that in both the cases (i) and (ii), Ai -I- A 2 > 0, we must have 
Ai -I- A 2 = A/|M| > 0. Here we assume A > 0 without loss of generality. Putting this result into Eq. (9), the 
dispersion of the surface states is finally given by 

Aarc(fcx, ftz) = sgn{M)Akx. (14) 


( 12 ) 

(13) 


The corresponding wavefunction can be simplified as 


(r) = 


sgn(M) 

1 


c>'iy 


— e 


>'2y\ 


(15) 


where (7 is a normalization factor and Ai _2 = A/2\M\ =F y^~(A/ 2 M)^Ar~^^ At the surface of y = 0, the surface 
states are eigenstates of ax with a uniform effective velocity, z;eff = sgn{M)A/ti. Thus they are chiral surface states 
around the surface parallel with the z direction. Also in both cases (i) and (ii), we have A 1 A 2 > 0 and henceforth 
Afc > 0. Therefore the solution of Eermi surface states is restricted inside a circle defined by k"^ + k^ < k^, as shown 
by Figs. 1(b) and (c). At zero Fermi energy, i.e., kx = 0, the surface states exist for all \kz\ < kc which produces a 
Fermi arc connecting two Weyl nodes. For a non-zero Fermi energy, the ends of the Fermi arc are shifted away from 
the Weyl nodes until they vanish. 
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FIG. 2: Top left: the energy spectrum of a Weyl semimetal as a function of at = ky = 0. Top right: the Landau bands 
of the Weyl semimetal in a z-direction magnetic field (B = 0.1 T). The 0th Landau bands are in red. Bottom: the DOS in the 
absence and presence of the z-direction magnetic field. The circles indicate the DOS at the Weyl nodes, which grows linearly 
with the magnetic field. Parameters: kc = 0.1/nm, M = 5 eV-nm^, and A — 1 eV-nm. 


III. FIELD-INDUCED DENSITY OF STATES AND CHIRAL ANOMALY 


After showing that the two-node model is capable of capturing all the topological properties, now we are ready 
to demonstrate its transport properties arising from its chiral properties. It is well known that the DOS of a Weyl 
semimetal vanishes at the Weyl nodes, following a E'^ dependence. This can also be captured by the two-node model 
in which the DOS at small energy E is given by 


N{E) = 


E^ 

2Tr^ A^Mkc' 


(16) 


as shown in Fig. 2 (see Appendix A for the general formulas of the DOS). When the Fermi energy locates right at the 
Weyl nodes, the DOS is exactly zero. At this exact semimetal phase, it is expected to have exotic electronic transports 
at zero temperature. It is found that a finite DOS can be generated by a magnetic field B along the direction that 
connects the two Weyl nodes, i.e., a z-direction magnetic field in the present work. The z-direction magnetic field can 
split the Dirac cones into a bundle of ID bands of Landau levels dispersing with kz [19, 40, 54], as shown in Fig. 2. 
Figure 2 shows how the DOS in the absence of the magnetic field evolves into a set of diverging peaks in the presence 
of the magnetic field. The divergence is due to the van Hove singularity at the band edges of the ID Landau bands. 

In the presence of a magnetic field in the z direction, the energies of electrons form a set of Landau bands dispersing 
with kz as following [40] 


Elf = uj/2 ± v>l 

Ef = uj/2- Mkl+ Mkl, zz = 0 (17) 

where Mi, = Mk^ — Mk1 — wz/, t] = 'J^Ajls and w = 2M/£^. Ib = \Jhl\eB\ is the magnetic length. Each Landau 
band has the Landau degeneracy = l/27r£^ = |ei?|//i in a unit area in the x-y plane. 

In the magnetic field, the DOS near the Weyl nodes is given by 

~ VM^yE + Mk^-uj/2' 


(18) 
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which is contributed by the 0th Landau bands. At the Weyl nodes, the DOS reduces to Nk^ = (l/27rMA:c) x (l/27r^^), 
which grows linearly with the Landau degeneracy Nl = \eB\/h and hence the magnetic field B. This resulting DOS 
leads to a semimetal to metal transition, which is under intensive investigation in the context of the disorder-induced 
DOS [55-59], instead of in the presence of magnetic field. 

In the z-direction field, the Fermi energy crosses with the 0th Landau bands at kc and —kc, where the Fermi 
velocities are opposite, i.e., vp = ±2Mkc/k. This chiral property of the 0th bands is exactly the scenario required 
in Nielsen and Ninomiya’s proposal for the chiral anomaly [1]. As an electric field is applied along the Weyl node 
direction, the changing rates of charge carriers are dN±/dt = ±.{evF near the two Weyl nodes at k^ = Eke, 
respectively. Thus charges can be pumped from near kc to —kc, literally leading to the non-conservation of chiral 
charge, i.e., the chiral anomaly [31]. Later, we will focus on the electronic transport in this situation when the chiral 
anomaly happens. In the following discussions, we will constrain to the case near the Weyl nodes, i.e., the Fermi 
energy is located between two Landau bands of 1±, and at very low temperatures, kpT <C w,? 7 . 


IV. LONGITUDINAL MAGNETOCONDUCTIVITY 


At sufficiently low temperatures, i.e., kpT uj, rj, and not far away from the Weyl nodes, the electronic transport 
can be effectively conducted by the 0th bands of Landau levels. When the electric and magnetic fields are in parallel 
with each other, the changing rate of density of charge carriers near one node is maximal according to the picture of 
chiral anomaly. In this case, the semiclassical conductivity of the 0th Landau bands can be found with the help of 
the standard Green function formulism [40] . Alternatively, it can be simply figured out by using the Einstein relation 
(jzz = e^NpD, where the DOS can be found as the Landau degeneracy times the DOS of one-dimensional systems, 
i.e., Np = (l/27r^^) x (l/Trhvp). D = is the diffusion coefficient in one dimension, is the transport time, 

Ui? is the Fermi velocity in the z direction and kp = hvp/2M is the Fermi wave number. For the scattering among 
the states on the Fermi surface of the 0th Landau bands, the transport time can be found as 


h 



K,K 




Vp 


(19) 


where j., represents the scattering matrix elements and (...) means the impurity average (see Appendix C 

for details). 

The transport time is sensitive to the scattering potential in materials. One of the convenient choices is the 
random Gaussian potential 


C7(r) 


E 




( 20 ) 


where Ui measures the scattering strength of a randomly distributed impurity at Ri, and d is a parameter that 
determines the range of the scattering potential. The Gaussian potential allows us to study the effect of the potential 
range in a controllable way, which we find it crucial in the present study. Now we have two characteristic lengths, 
the potential range d and the magnetic length ip, which define two regimes, the long-range potential regime d^ ip 
and the short-range potential limit d ip. Note that, for a given d in realistic materials, varying the magnetic field 
alone can cross between the two regimes. Empirically, the magnetic length ip = 25.6 nm /^/\B\ with B in Tesla. 
In the strong-field limit, e.g., R > 10 T, the magnetic length ip becomes less than 10 nm, it is reasonable to regard 
smooth fluctuations in materials as long-range. 

With the random Gaussian potential, we can find the transport time as well as the conductivity. In particular, at 
the Weyl nodes the transport time is obtained as (see Appendix B 1 for details) 






imp 


.0’“- 2TTMkc 2d2 -p £2 


( 21 ) 


and hence the longitudinal conductivity 


UzziB) 


{2Mkcf{2d? 




V i^ 

^ imp^g 


BJ 


( 22 ) 


where Idmp = measures the strength of the scattering and V = L^LyL^ is the volume of the system. L^^y^z 

are the sizes of the system along the x, y and z directions, respectively. This conductivity is generated by the inter¬ 
node scattering with a momentum transfer of 2kc. As the magnetic field goes to zero, the magnetic length diverges 
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FIG. 3: The longitudinal conductivity a^z and transverse conductivity of the Weyl semimetal in the z-direction magnetic 
field B for different potential ranges. The shared parameters: kc = 0.1/nm, M = 5 eV-nm^, A = 1 eV-nm, Fimp = 10 
(eV)^-nm®. 


and d/is 0, and Eq. (22) gives a minimum conductivity 


CTzz(O) 


h Viinp 


(23) 


even though the DOS vanishes at the Weyl nodes at zero magnetic field. A similar result was found in the absence of 
the Landau levels [58]. 

According to d, we have two cases. (1) In the short-range limit, d = 0, then a^z does not depend on the magnetic 
field, giving a zero magnetoconductivity, which recovers the result for the delta potential [40, 41]. (2) As long as the 
potential range is finite, i.e., d > 0, we can have a magnetoconductivity. Using Eq. (22), 


Aa^^{B) 


<Jzz{B) - CTzziO) _ 

CTzz(O) Bo’ 


(24) 


where Bq = H/2ed^. Thus the magneto conductivity is given by the range of impurity potential, and independent 
of the model parameters. This means that we have a positive linear z-direction magnetoconductivity for the Weyl 
semimetal. A finite carrier density uq can drive the system away from the Weyl nodes, then kc in Eq. (22) is to be 
replaced by kp = kc + sgn(M)27r^£^no. The finite Uq can vary the linear-B dependence, but a strong magnetic field 
can always squeeze the Fermi energy to kc, and recover the linear magnetoconductivity. 
















7 


A linear-i? magnetoconductivity arising from the Landau degeneracy has been obtained before [37, 39], based on 
the assumption that the transport time and Fermi velocity are constant. However, in the present case, we have taken 
into account the magnetic field dependence of the transport time, and thus the B-lmeai magnetoconductivity here 
has a different mechanism as a result of the interplay of the Landau degeneracy and impurity scattering. Also, in 
the presence of the charged impurities, a magnetoconductivity can be found in the quantum limit [41]. A B^ 
magneto conductivity can also be found in the semiclassical limit [36, 37]. 


V. TRANSVERSE MAGNETOCONDUCTIVITY 


When electric and magnetic fields are perpendicular to each other, the changing rate of density of charge carriers 
near each node vanishes. In this case, because the Landau bands in the z-direction magnetic field only disperse with 
kz, the effective velocity along the x direction Vx = dEo/hdkx = 0. The leading-order x-direction conductivity arises 
from the inter-band velocity and the scatterings between the 0th bands with the bands of 1±, which are higher-order 
perturbation processes. Thus the transverse conductivity is usually much smaller than the longitudinal conductivity. 

When the Fermi energy crosses only the 0th bands, the leading-order conductivity is given by [60] 

cTxxiB) = a+^{B)+ a-^{B), 

^ E ^«kyo,i±Gt±,kyi±,o) , (25) 

where = ^/{Ep — E^^ + is the retarded Green’s function of electrons in the 0th bands, = 

l/{Ep =F EI"^ — ih/2Tl^) are the advanced Green functions of electrons in the 1± bands, respectively, and 
are the corresponding lifetimes. The inter-band velocities along the x direction Vg are given by 


V 


X 

0 , 1 + 


V 


X 

0 , 1 - 


(-B ( 


COS ■ 




■ uj sin • 


y 


y/2h\ 


r] sin 


y_ y 

2 2 


(26) 

(27) 


where cos 9l = Mij \/M.\ + rj^. Considering that Fermi energy crosses only the 0th bands and that h/r^ and k/r]^ 
are very small, Eq. (25) can be simplified to 


yiB) = 


~V~ 


E i<i± 


S{Ep-El) h 
[Ep - Ely 2rlf ■ 


(28) 


The lifetimes due to the random Gaussian potential can be calculated by using the Born approximation, 

h 


A± 


K,K 


After lengthy but straightforward calculations, they can be obtained as (see Appendix B 2 for details) 


h 1 + 


(1Tcos6»7). 


Wnhvp {(P + G^/2)'^ 

The substitutions of Eqs. (26) and (27) for Uq and of the results (30) for into Eq. (28) give 


<Jxx{B) - {2d^ + iy 


where 


Hkz)= ^ 




yi-Eyf ■ yi-yy 




(29) 


(30) 


(31) 


(32) 












When the Fermi energy is at the Weyl node, the transverse conductivity Uxx at the Weyl nodes is given by 




(? Mmp £%(l + e ^ 

h WTT^M^k^ (2d2+£|)2 


( 33 ) 


where J' = (2uj^ + 7}^)/(4oj^ + 4r/^), so « 1/2 in the strong-field limit (£b M/A) and « 1/4 in the weak-field 
limit {£b ^ M/A). We choose the value of l/mp so that the band broadening is always much weaker than the spacing 
between the 0th and 1st Landau bands for an arbitrary magnetic field. It is safe to assume that only the 0th Landau 
band contributes to the conductivity at half filling. There are three cases as shown in Figs. 3 (d)-(f). At d = 0, axx 
reduces to the result for the delta potential and tr^x oc S, a linear magnetoconductivity as Uzzi but much smaller [40]. 
In the long-range potential limit £b, we have cr^x ~ l/^i which gives a negative magnetoconductivity. For a 
finite potential range d, we would have a crossover of Uxx from i?-linear to l/B dependence. Alternatively, as shown 
in Fig. 3 (e), for a finite d (= 5 nm) comparable to the magnetic length £b^ we have a crossover of Uxx from a linear-i? 
dependence in weak fields to a 1/B dependence in strong fields. While at d = 0 and d^ £b, we have the two limits 
as shown in Figs. 3 (d) and (f), respectively. For shorter d, a larger critical magnetic field for the crossover is needed. 
Figure 3 also shows that the conductivity is larger for shorter d, so the 1/B transverse magnetoconductivity in the 
long-range limit may not survive when there are additional short-range scatters. 

In particular, in Fig. 3 (f), axx oc 1/B in the long-range potential limit. In the field perpendicular to the x-y 
plane, there is also a Hall conductivity ayx = sgn(M)(fcc/7r)e^//i -I- eno/B, where the first term is the anomalous 
Hall conductivity and the second term is the classical conductivity. In weak fields, the classical Hall effect dominates, 
then both Gxx and Oyx are proportional to I/H, and the resistivity pxx = o’xx/(ct^x + ^yx) is found to be linear in B. 
Note that here the linear MR in perpendicular fields has a different scenario compared to the previous works [61, 62]. 
Abrikosov used the Hamiltonian vk ■ a with linear dispersion and modelled the disorder by the screened Coulomb 
potential under the random phase approximation [61]. Song et al. discussed a semiclassical mechanism [62]. 


VI. PARAMAGNETIC WEYL SEMIMETAL 

The existence of the non-zero Chern number or chiral surface states indicates the time reversal symmetry breaking 
in the two-node model in Eq. (1). Correspondingly a quantum anomalous Hall conductance appears. To have a 
paramagnetic Weyl semimetal, we have to introduce a time reversal counterpart for the two Weyl nodes in Eq. (1). 
A straightforward extension is as follows 

H = A{kxax + kyay) + M{k/- k‘^)P, (34) 

where the Dirac matrices are oix = ax ® ax, oty = ax ® ay^ j3 = a^ ® az- It contains four Weyl nodes, which are 
doubly degenerate. The surface electrons around the z direction consist of two branches with opposite spins and 
opposite effective velocities. It will give rise to the quantum spin Hall effect, compared to the quantum anomalous 
Hall effect in a Weyl semimetal of a single pair of nodes. The dispersions of two branches of the 0th bands are 
E±{kz) = ±(w/2 — Mk/ + Mkl). In a weak magnetic field, the magnetoconductivity is also linear in magnetic field, 
as that for Weyl semimetal in Eq. (1). However, for a strong field, the strong field dependent Eermi energy will give a 
different field dependence of conductivity as the two 0th bands shift away from each other as increasing field. Einally 
it is worthy pointing out that this kind of semimetals is different from a simple Dirac semimetal described by a pair 
of degenerate Dirac cones, Evk ■ a (e.g., at the phase transition point between a 3D topological insulator and trivial 
insulator), in which there does not exist the surface states as the two nodes are not separated in momentum space 
and the model is topologically marginal. Also, the Zeeman effect may also contribute to a linear contribution to 
the magnetoconductivity [63] . But its sign depends on the g-factor of sample and usually its magnetoconductivity is 
negative, opposite to that in this work. 


VII. CONCLUSIONS AND DISCUSSION 

The key conclusion of the present work is the positive and linear magnetoconductivity of Weyl semimetals near the 
Weyl nodes. The Fermi energy is assumed to cross only the 0th Landau bands of the semimetal at low temperatures 
{Ep u) and fcsT <C cu). The magnetoconductivity depends on the types of scattering potentials, in which the 
potential range is a characteristic parameter. Our conclusion is different from the theoretical predictions of the positive 
B^ magnetoconductivity in several previous works [36, 37], because where higher Fermi energies (Ep ^ kpT) and 
^ ^ were assumed, therefore there were many more Landau bands on the Fermi surface. Recently, a positive 
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longitudinal magnetoconductivity has been observed in several different candidates of topological semimetal [44-51], 
and was claimed to be related to the chiral anomaly. We believe that the linear term is one of the indispensable 
ingredients in the formula of the magnetoconductivity in the intrinsic Weyl semimetal phase. When more Landau 
levels come into play as the Fermi energy is shifted from near the Weyl nodes, the magnetoconductivity is expected 
to deviate from linear. 

The result of the magnetoconductivity is based on the Born approximation. When the magnetic length becomes 
much shorter than the range of the disorder potential, electrons may be scattered by the same impurity for multiple 
times. The Born approximation contains the correlation of two scattering events by the same impurity [64]. In this 
situation, the validity of the Born approximation was questioned in two dimensions [65, 66]. In three dimensions, it is 
still unclear whether the correlation of two scattering events in the Born approximation is the building block for the 
multiple scattering under extremely strong magnetic fields [62, 67]. So far, theories in three dimensions employ the 
Born approximation, e.g., the quantum linear magnetoresistance [61]. The treatment beyond the Born approximation 
will be a challenging topic for three-dimensional systems under extremely strong magnetic fields. 
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Appendix A: Density of states 


1. In the absence of magnetic field 


In the absence of magnetic field, the DOS can be calculated as N{E) = ~ After the 

straightforward calculations, the results of DOS are summarized as followings. 

If > 2M^fc^, the results of DOS are given by 


^ \v 2{E), \E\ > Mkl 


(Al) 


where 

^ ^ 2A^M{Mkl+E)+A^ + 2ME 

2A^M{Mkl -E)+A'^- 2ME 
2A^M{Mkl + A) + + 2ME 

^ ^ y/AM^E^ +A^- AA^M^kl 

If A^ < 2M^fcg, the DOS is instead given by 

^ |£;|<AV2M 

[d2(A)+D4(A), \E\>Mkl 

where 

, , , AJAM^kl - 2A2 -h v'AM^E^ - A* 

T); 2 {E) — In- , =-, 

V'4A2M2fc2 _ - AM^E^ 

, , , AJAM'^kl - 2A2 -t ^AM^E^ - A^ 

VAE) = In—— ■ r-. 

^/AM^E^ -h A4 - AA^M^kl 


(A2) 

(A3) 

(A4) 

(A5) 

(A6) 


For a small energy E ^ Mkl, Al' I2M (in both the cases A? > 2M^k^ and A^ < 2M^fc^), the DOS can be simplified 
as 


N{E) = 


27r2A2MA:c‘ 


(A7) 



















10 


2. In the presence of magnetic field 


In the presence of magnetic field, the DOS can be calculated as N{E) = (l/27rfgLz) s=± ~ ^)- Sum¬ 

ming over kz and s = ±, N(E) can be written as 


N{E) = 


\E-uj/2\ 


N. 

E 


+ 1 


% V^E) V Mkl -iyLo + V,{E) 


N. 


+ E 


1 


1 


.=0 V,iE) V Mkl - I/W - V,{E) 


(AS) 


where Vu{E) = ^/{E — w/2)2 — i/jf' and N± are the largest integers for which Mk'^ — vu) ± Vv{E) are positive, 
respectively. 

In the quantum limit, only the v = 0 term is retained and N{E) reduces to 


47r2£| y/M^E + Mkl-uj/2' 
which is approximately proportional to the magnetic field. 


(A9) 


Appendix B: Lifetimes and transport times 


1. Transport time r. 


0,tr 


The transport time at the Fermi surface is calculated as 


h 

tr 


7 / 7 / 


(Bl) 


KK 


Using 5{Ef — E ^,) = [S{kF — k'^) -I- 6{kF + k'^)]/hvF and substituting the expressions Eq. (C16) for the scattering 
matrix elements p), 


fc'.fc'E ^ ^ 


hvF 


x5{q, +ki- k,)[S{k', + kF) + 5{k', - fcF)](I - ^). 

Vp 

where q\ = q^ + qy. Changing the summations over k'^, fcf, to integrals, Eq. (B2) can be simplified as 




The substitution of the Fourier transform of the potential Ui{q) = UiC ? /2 gives 

^ _ ^^imp I" 


n. 


hvF J (27r)2 


where Ump = Ei /^- Performing the integral in polar coordinates and we finally obtain 


Ui. 


^-4d^k% 


hvF 2 f{iP E ^b/2) 


(B2) 


(B3) 


(B4) 


(B5) 
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2. Lifetimes r, 


i± 


We assume that the impurity scattering is so weak that when the Fermi energy crosses only the 0th landau bands, 
the leading order of the lifetimes of the states in the 1± Landau bands can be found as 


-I± = !*)«(£,-E“,). 


(Be) 


k' ,k' 


Substituting Eq. (C23) for the scattering matrix elements ^ P) and then using 8{Ep — E °,) = [S{kF — 

k'^) + 5{kp + k'^)\/hvp, we have 


h 

‘kp 




kx ^kfr^k'^ ,k'^ 

dqxdqy 9 /2 

2 


-00 (27r) 


^ [(^ {Qx ^ qy j imp E {U (^qxT qy^^^kp'j') 


impj 


(B7) 


Substituting the Fourier transform of the Gaussian potential and then performing the integrals, Eq. (B7) become 


h 

3± 


H 


WiThvp {d'^+i'^/2)^ 
3e 

Vimp^B 1 + 


(1TCOS0^^). 


(B8) 


In the strong-held limit £b —t 0, cos 9l^ —>■ —sgn(M). Suppose M > 0, then in the strong-held limit, h/rj.^ = 0 while 


Snhvp {d?+£\/2Y 


(B9) 


Appendix C: Scattering matrix elements 


To evaluate the square of the scattering matrix elements under the average of impurity conhgurations, we need to 
calculate The sum runs over all impurities. The integrals is dehned as 




L^L: 


I dr^i:{y)ip^,,{y)U,{r - R^)E(K-k.)x+^(k'^-k. 


(Cl) 


where t/i(r —R^) denotes the scattering potential of a single impurity at the position Rj. The y-direction wavefunction 
^kSy) given in Eq. (Cll). Fourier transforming the potential by using 


C/,(r-R,) = 




(C2) 


and then using the formula dxe^^^ = 2pS{k), can be rewritten as 


7,y,p(Rj) = 


d^ci ^-»q.R, 
Stt LxL z 


'’C/,(q) / dyLpll{y)ip'^, {y) 


- kx)5{q^ +K- ^z)- 

Using Eq. (C3), can be written as 


(C3) 


'i-J 

= X] / 2nL\ ^~"'^'^"Ui{ci)S{qx + fc(, - kx)S{qx + k'^ - k^) 

-j J X Z 

X j dyipll{y)(p'^,^iy)E‘‘^y 

X I ^^^^e^^'-^^UM')S{q'x + K-kx)S{q', + k',-k.) 

X J dy'ipiSy'M;{y')e-<y'. (C4) 
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Under the average of impurity configurations, we employ the theory of impurity average [64] 



‘•^*C7.(q)C/,(q')) « 5 ] t/^(q) Vq' = (27r)3(C/2(q))i„pd(q - q'), 

i 

(C5) 

where {U^{(l))imp = Ei Uf{(i)/V. 

Then we can integrate over q' to obtain 


7 - 1 /',^' 




= j 2-kL L + kx)5{qz + k'^ kz) 



X J dyipi:{y)yl,^{y)e^‘^yy j dy'yt{y')yi,* {y')e-^‘^yy', 

(C 6 ) 

where we have used 

d'^{qx + k'x-kx) = ^S{qx + k'x-kx), 

(C7) 


S'^iqz + kx - kz) = ^6{qz + k'z - kz). 

Ztt 

(C 8 ) 


1. Scattering matrix element ,y y p) 


In the Landau gauge A = —Byk, the corresponding eigen wavefunctions for v > \ can be written 


while for = 0 as 


{v\v+,kx,k:,) = 
{T\v-,ka:,kz) = 

{v\Q,k:^,k^) = 


QV 

cos 

sin-^ 

9^ 

sin 

9^ 

-cos- 5 ^ 1pu,k^kA^) 


0 

i’o,k^kA^) 


where cos 9'^ = M.^!\/M-t + v-rf. The wavefunctions '4’v,k^kA^) given by 

^ikxX-\-ikzZ 


'4^v,kxkz (l*) — 


^y,kAy) = 


V LxLz 


-'yv,kAy)^ 






(C9) 


(CIO) 


(Cll) 


where dq = kxi% is the guiding center and T-Li, are the Hermite polynomials. 

Using the wavefunctions Eq. (C9), the scattering matrix elements between state |0, kx, kz) and state |0, fc(,) can 

be written as 


uZkyk'k' ^ (0,fc.,fc,|[/(r)|0,fc;,fci)=E,4.o(R.)- 


and thus .y y p) = Tg q . With v = 0 and /r = 0 in Eq. (C 6 ), we have 


.7-0.0 

-^ 0.0 




QiTtL^L/z 


{U'^{q))inipS{qx +k'x- kx) 


x5{qz + k'z- kz) / dyipll {y)(pl^-q^ 


(C12) 


(C13) 
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Substituting the explicit form of the wavefunction Eq. (Cll), we can find 


and hence 

Therefore Eq. (C13) becomes 

{\Ul\ ■.'.'?) = 


dy^kAy)^k^-g^ 


J 2'kL L + k'^ — kx)5{qz + — kz). 




(CM) 


(C15) 


(C16) 


2. Scattering matrix element .y y p) 


Similarly, the scattering matrix elements between state |1±, kx-, kz) and state |0, fc(,, fcf,) can be written as 


0l 

sin^/i,o(R*), 

(C17) 

Ql 

cos ^/yo(Ri), 

(CIS) 


and thus y p) = T cos0^^)/2. With v = \ and /r = 0 in Eq. (C6), we have 


.7-1,0 _ 
-^1,0 ~ 


/rr2 


{U (q))imp / dyipf,l{y)(pl, (y)e^iyy 


J QiTtLx^z j 

^x)^{,Qz “1“ ^z)- 


Substituting the wavefunction Eq. (Cll), we can find 

/ dy(pl*(y)ipl_„(y)e^‘>^y = 


_ (-Bjqx + 'i'dv) -e%(q^, -2iq^q^+iAk^q^)/A 

V2 


(C19) 


(C20) 


and hence 


dy<Pk*Aym,-Ay)e"‘^'‘^ 


?i^p-9i4/2 


Therefore Eq. (C19) can be rewritten as 


p2 

T-1,0 _ ^ 

^ 1,0 - 2 


J 2'kL L {A))d{qx + fc(, — kx)S{qz + — kz), 


and 




xS{qx + k'^- kx)5(qz + k'^- kz). 


(C21) 


(C22) 


(C23) 
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